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Abstract. We compute the triple gluon, quark-gluon and ghost-gluon vertices of QCD at 
the symmetric subtraction point at two loops in the MS scheme. In addition we renormalize 
each of the three vertices in their respective momentum subtraction schemes, MOMggg, MOMq 
and MOMh. The conversion functions of all the wave functions, coupling constant and gauge 
parameter renormalization constants of each of the schemes relative to MS are determined 
analytically. These are then used to derive the three loop anomalous dimensions of the gluon, 
quark, Faddeev-Popov ghost and gauge parameter as well as the /3-function in an arbitrary 
linear covariant gauge for each MOM scheme. There is good agreement of the latter with earlier 
Landau gauge numerical estimates of Chetyrkin and Seidensticker. 
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1 Introduction. 



The structure of the renormalization group functions of Quantum Chromodynamics (QCD) has 
been established to a high degree of precision for well over a decade now. Originally the one 
loop discovery of asymptotic freedom due to the negative (3- function, [HE], was swiftly followed 
by the two loop computation, [3l Hj. Within a decade the three loop term emerged, [5], but 
the four loop result, [6], took substantially longer to determine. It was later confirmed in [7]. 
To appreciate how involved the calculation of [6] was, it required the evaluation of the order of 
50000 of Feynman diagrams and an intense amount of symbolic algebraic manipulation using 
the language Form, [8]. Also new techniques to evaluate four loop Feynman integrals were 
developed, [8], as well as the coding of a Form routine to handle the colour group algebra 
manipulations automatically, [3]. In addition to the /3-function governing the running coupling 
constant the other renormalization group functions have been available at a variety of loop 
orders over the same time scale, [10} ITT| \T7[ \13\ [l~4"t I15j. In summarizing the status of the basic 
renormalization group functions of QCD we have concentrated on the status of the modified 
minimal subtraction, MS, scheme. This renormalization scheme is the most commonly used 
for perturbative QCD computations. It is used primarily because it is a mass independent 
scheme. Consequently one simplifying feature is that all the Feynman graphs which one needs 
to evaluate to determine the /3-function to four loops in this scheme involve single scale 2-point 
Feynman integrals. This means, for instance, that techniques such as integration by parts can 
be systematically used to reduce integrals to either basic master topologies or simple chain 
type integrals whose evaluation is effectively trivial. Indeed the Mincer algorithm, [16], is an 
excellent example of the implementation of integration by parts and has been coded in Form, 
[8]. Briefly Mincer determines massless 2-point integrals at three loops to the finite part with 
respect to the regulator. The divergences are written as simple poles in e where d = 4 — 2e 
and d is the arbitrary spacetime dimension of dimensional regularization. The latter, which 
we will use, is the main regularization method for perturbative quantum field theories and, 
for instance, preserves gauge invariance. Whilst Mincer has been used extensively for many 
problems other than the basic renormalization group functions of QCD, the FORM version, |17| . 
was used to renormalize the QCD coupling constant at three loops, [12J. The reason for this 
is that one can nullify an external momentum of one of the legs of a 3-point function, without 
introducing spurious infrared infinities, and hence apply the Mincer algorithm to obtain the 
correct three loop MS /3-function of [5j. Although this has proved to be a powerful technique 
for renormalization group functions, for more physical problems the MS scheme has several 
drawbacks. 

One of these is the fact that the MS scheme is not a physical scheme. Specifically, the 
subtraction at the vertex as noted earlier in the Mincer approach is at a point of exceptional 
momentum. Whilst this is sufficient to extract the divergences and hence find the renormal- 
ization constant, one would have potential infrared ambiguities if one were trying to make a 
measurement of the finite part of the associated vertex function at this exceptional momentum 
configuration. This has been known for a long time, [18} 119] . but has also been remarked upon 
again more recently in the specific context of Green's functions with an inserted operator, [20] . 
Indeed the measurement of vertices is a topic of interest as one can study any of the vertices 
non-perturbatively using say lattice regularization. Such analyses are necessary for determining 
the strength of the strong coupling constant accurately for comparison with experiment and 
examining its behaviour at low energy. Therefore, it seems appropriate to consider vertex mo- 
mentum configurations which are non-exceptional and hence schemes which are physical rather 
than the unphysical minimal schemes. A set of such schemes was introduced in [21] and given 
the general designation of momentum subtraction schemes being denoted by MOM. In essence 
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they consist of two criteria. The first is that 3-point vertices are considered at a symmetric 
subtraction point, [21j . In other words the squares of the three external momenta are all set 
equal to each other. Therefore, there is no nullification of an external momentum and hence 
no exceptional momenta. The second is that at the subtraction point the scheme is defined in 
such a way that after renormalization there are no 0(a) corrections where a = g 2 /(16ir' 2 ) and 
g is the gauge coupling constant. Thus the renormalization constants all contain finite parts, 
|21j . in addition to the poles in e which must always be subtracted. In [21] the MOM schemes 
were analysed comprehensively at one loop. By schemes we mean the three types derived from 
the respective three 3-point vertices in the canonical linear covariant gauge fixing. These are 
denoted by MOMggg, MOMh and MOMq due to their origin from the respective triple gluon, 
ghost-gluon and quark-gluon vertices. Indeed these one loop MOM scheme computations of the 
renormalization group functions have been the state of the art for a long time for a relatively 
simple reason. This is because for the vertex renormalization the symmetric subtraction point 
momentum configuration introduces single scale 3-point Feynman integrals, [21J. At one loop 
there is only one such basic master integral to evaluate which was given in [21]. However, at two 
loops there are several basic master integrals which were only evaluated analytically in recent 
years. Therefore, it is the purpose of this article to extend the work of [21 J to two loops. We will 
achieve this by evaluating all 2-point and 3-point vertices to the finite part at two loops at the 
symmetric subtraction point analytically. Hence we will determine the gluon, Faddeev-Popov 
ghost and quark wave function and coupling constant renormalization constants for each of the 
three MOM schemes in an arbitrary linear covariant gauge. En route we will also provide the 
structure of each of the three vertices to the finite part in the MS scheme at the symmetric point. 
This information should prove useful to lattice groups seeking to measure any of the vertex func- 
tions in the MS scheme. For instance, in such an exercise the lattice results must match onto 
the ultraviolet part of the vertex function which is where perturbation theory is valid. Equipped 
with the two loop renormalization constants we will also determine the conversion functions for 
each scheme including the relation between the coupling constants. Hence via properties of the 
renormalization group we will determine each MOM scheme renormalization group function to 
three loops in an arbitrary linear covariant gauge. 

Finally, we mention other related work in this area. Prior to this article there were approxi- 
mate calculations of the three MOM /3- functions in the Landau gauge, [22]. This was achieved 
by approximating the basic master two loop integrals at the symmetric point by an expansion 
where one of the external momentum is marked to produce an asymmetric Green's function. 
Then the integrals were expanded in the ratio of this marked momentum to the other indepen- 
dent momentum, |22[ 123] . So if a sufficient number of terms is computed in this parameter, 
where the coefficients are rapidly decreasing in size, then a reasonable approximation can be 
deduced by truncating at an appropriate order, [22], [23] . Moreover, error estimates can be de- 
duced. In this expansion the 3-point functions reduce to 2-point functions. Hence each term 
in the series is evaluated using the Mincer algorithm. Indeed we previewed the results of this 
article in [24J, by providing the exact three loop QCD MOM /3-functions in the Landau gauge 
and demonstrated how accurate the results were in comparison to |22j . This was a very close 
overlap which was impressive given the level of computing technology available to the authors of 
[22j at that time. Next, it would be remiss if we failed to mention a related three loop MOM (3- 
function computation given in [25]. There the MOM /3- function was defined using the invariant 
charge concept of [261 127j . As it involves only finite parts of 2-point functions it is independent 
of any of the vertices of the theory unlike the MOMggg, MOMh or MOMq schemes. Hence, 
aside from the Riemann zeta series, it does not have any of the special number structures, such 
as harmonic polylogarithms, which derive from the single scale 3-point master integrals and 
are evident in our analytic results, |24j . Also in this context a variation on this theme of es- 
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tablishing MOM scheme definitions based solely on 2-point functions has been developed more 
recently in [28]. Known as the minimal MOM scheme it circumvents a renormalization condition 
based on the ghost-gluon vertex 3-point function by imposing the alternative condition that the 
renormalization constant of the ghost-gluon vertex is the same as that of the vertex in the MS 
scheme itself. One benefit of this is that in a linear covariant gauge a non-perturbative run- 
ning coupling constant can be defined purely in terms of the gluon and ghost propagator form 
factors. With this definition the running coupling constant can be measured more accurately 
in principle using lattice gauge theory techniques. For instance, recent activity on this specific 
aspect can be found in [29] [30]. Given that one can determine the gluon and ghost 2-point 
functions to several loop orders in perturbation theory, the four loop minimal MOM /3-function 
has been determined in [28]. Though unlike [25J quark mass effects have only been estimated 
in |28j . Finally, we note that the measurement of vertices non-perturbatively is not exclusively 
studied by lattice techniques. For instance, recently the triple gluon vertex was examined using 
Schwinger-Dyson methods, [31] • Therefore, the structure of the two loop vertices given at the 
symmetric point here should also prove relevant in analyses using those techniques too. Other 
two loop studies of 3-point vertices of QCD include the results of [32], [33] where the triple gluon 
vertex was examined in the zero momentum limit and in the on-shell configuration respectively. 
The latter also includes the same analysis for the ghost-gluon vertex. 

The article is organized as follows. We outline the general formalism and notation we use in 
section 2 as well as discussing various aspects of the renormalization in each of the three MOM 
schemes in section 3. The results for the MS and MOM amplitudes of the respective vertex 
functions as well as the conversion functions and the three loop renormalization group functions 
are given in each of the three following sections^. Finally, we present conclusions in section 7. 
An appendix contains the explicit tensors of the bases for each of the three vertices as well as 
the respective projection matrices. 



2 Preliminaries. 

We begin by discussing the general features of the computation we perform. The three Green's 
functions we consider are (^A«(p)Ab(q)A%(r)\, {^{pW \q)A%{r)) and (c a (p)c b (q)A%(r)^ where 
r = — p — q by momentum conservation. The symmetric subtraction point is defined by the 
condition 

2 2 2 2 /„ 

p = q = r = — [i (2.1) 

where fi is the common mass scale. It will also be used as the mass scale to ensure that the 
coupling constant remains dimensionless in dimensional regularization in d-dimensions which we 
use throughout. Therefore our results for the finite parts of the vertex functions will not involve 
logarithms which can be restored from knowledge of the renormalization group functions. From 
(12. ID we have 

pq = \f (2.2) 

and we will use p and q as the two independent momenta throughout. Their sum will be taken 
to flow out through a gluon external leg which will therefore be the reference leg in each of the 
vertices. Given that each vertex has a colour group tensor associated with it, we factor it off 
when we consider the symmetric point, which we do exclusively from now on, by defining 

A«(p)AUq)A c a (-p - q)) 2 _ 2 __ 2 = f abc Z®%(P,< 



*A11 the results presented in the article, including the full analytic forms for an arbitrary gauge, have been 
included in an attached electronic data file for each of the three MOM schemes. 
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^{p)f>(q)A c a {-p-q) 
c a (p)c b (q)Al(-p-q) 



(2.3) 



We will use ggg, qqg and ccg in equations to denote the triple gluon, quark-gluon and ghost-gluon 
vertex functions respectively. Next we decompose the Lorentz amplitudes £^fi§-(p, q)\ p 2 =q 2 = _ fl 2, 
Sa qs (p, q)\„2- n 2-_,,2 and T,a Cg (p, q)\„2-„2--,,2 into the scalar amplitudes. The Lorentz tensor 



basis for each function is not the same. Nor is the choice of basis we will use for the decomposition 
unique. The explicit forms of the tensors are given in Appendix A. Though we note that away 
from the symmetric point, where the equalities of (|2.ip are no longer valid, then the basis will 
involve a larger number of tensors. Therefore, we formally define the scalar amplitudes as 
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k=l 



k=l 
2 



(2.4) 



fc=i 



where £^(p, q) are the scalar amplitudes at (|2.1j) and V)^ 1 " iln% (p, q) are the tensors of the 
respective bases with i corresponding to one of ggg, qqg or ccg. We use k to label the basis 
elements and have chosen the labelling in such a way that channel 1 corresponds to the tensor 
which occurs in the Feynman rule of the corresponding vertex in the QCD Lagrangian. Though 
given the structure of the triple gluon vertex the first six tensors are part of that vertex. 

To determine the values of each of the scalar amplitudes we use the method of projection. 
In other words an identified amplitude can be isolated by multiplying the Green's function by a 
specific linear combination of the basis tensors. Therefore, we have 



r 6c Sff(p,g) 



\k) 
-qqg 

-ccg 

\k) 



Z5E«F(p,ff) 



M%* (V^ Sflva (p,q) (A-{p)Ai(q)A%(-p - q) 
(0 



M ccg ( v cc g 
Jvl kl {' (I) 



M g >, q) {^{pW{q)A c a {-p - q) 



(p,q)(c a (p)c b (q)Al{-p-q) 



(2.5) 



for each vertex where we have included the passive colour factor on the left hand side to comple- 
ment the one which is implicit on the right side. The free spinor indices in the quark-gluon vertex 
have been left implicit. The matrix A^L is the projection matrix and the explicit forms for each 
of the vertices are given in the appendix. It is computed by first finding the matrix J\fh for each 
vertex which is constructed from the basis tensors by Lorentz contraction in (i-dimensions using 
the conditions of the symmetric point, (I2.ip . In other words 



Ml 



k! 



(v. 



(k) fj,i...fj, n _ 



(p, q)V, 



i /ti.../t„ 

(0 



1 (p, q) 



(2.6) 



This produces a matrix, Mh, whose entries are polynomials in d and .ML corresponds to its 
inverse. For the quark-gluon vertex the tensor basis necessarily has to be built from 7-matrices 
in addition to the momentum vectors. As we will be working in dimensional regularization we 
will use the generalized 7-matrices, r^C"^ n , [341 1351 [36] . These form a complete set of matrices 
which span the spinor space of the associated d-dimensional spacetime. They are defined to be 
completely antisymmetric in the Lorentz indices and are given by 



1 (n) 



7 



[Mi 



■ r 



(2.7) 
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where an overall factor of 1/n! is understood and IVo) is the unit element. One beneficial 
property, among other general properties |37l 138] . is that the trace over the generalized 7- 
matrices is isotropic as 

tr (rg l y Mm r^- l/ ") oc s mn i^-^-^ (2.8) 

and l^i—ltmvi—vn \ s the unit tensor. For the quark-gluon vertex only the n = 1 and 3 generalized 
matrices arise, as we are working in the chiral limit throughout, which can be seen in the explicit 
decomposition in appendix A. 

We have used several main working tools to complete our analysis which are all computer 
based as it would be virtually impossible to proceed without automatic Feynman diagram gen- 
erators as well as symbolic manipulation programmes. For each of the three vertex functions 
the Feynman graphs are constructed with the Qgraf package, |39j . For the triple gluon vertex 
there are 8 one loop and 106 two loop diagrams contributing to the Green's function. For the 
other two vertices the number of graphs is the same with 2 one loop and 33 two loop diagrams. 
From the Qgraf output the Lorentz and colour indices are appended in the symbolic manipu- 
lation language Form, [8j. Indeed Form is used as the machinery for the rest of our algebraic 
computations as it is efficient in handling the huge amounts of algebra required for our analysis. 
To compute the Feynman graphs to the required finite part in dimensional regularization we use 
the Laporta algorithm, [3D]. Briefly the aim is to write each Feynman diagram in the Green's 
function in terms of a set of basic scalar master integrals whose expressions to the finite part 
are known. After applying the projection matrix to the Green's function the resulting scalar 
Feynman integrals are written in a specific format. Given the structure of the QCD propagators 
and vertices the scalar products in the numerators are rewritten in terms of the propagator 
denominators in preparation for using the method of [40J. However, given the symmetric point 
condition then for all the topologies there will be irreducible numerators where there are no cor- 
responding denominators. Equally there will be propagators raised to a power larger than unity. 
To reduce this very large number of scalar integrals to the set of masters requires an intense 
amount of integration by parts. One method which achieves this is the Laporta approach, |40| . 
This systematically determines all the integration by parts relations, as well as Lorentz identity 
relations, between all the integrals which are needed. The algorithm then uses a systematic way 
of reducing integrals classified in various levels to the lower levels or to a basic master in that 
level. The beauty of the method is that it terminates and can be coded for implementation on 
a computer. There are several available packages. We have chosen to use Reduze, [H], which 
is written in the symbolic manipulation formalism of GiNaC, [42j . whose working language is 
C++. One main aspect of the package which we exploit is to construct a database of relations 
between the integrals and then to lift out those we require for our specific computation. These 
are simply mapped to Form format and an integration module included within our overall au- 
tomatic Form programme. This sums up all the contributions to each of the Green's functions 
allowing us to perform the renormalization in any of the schemes of interest. For the latter 
we follow the algorithm for automatic Feynman diagram computations devised in [12]. This 
involves performing all the integrals as a function of bare parameters with the renormalized 
values introduced by a simple rescaling via the respective renormalization constants. Therefore, 
this means that we only compile the results for the vertex amplitudes once prior to determining 
the MS or MOM scheme renormalization group functions and amplitudes. 

Finally, as the analytic results we derive involve cumbersome expressions even in the Landau 
gauge as will be evident, we will give numerical values for all our results in an arbitrary gauge. As 
our computations revolve around the symmetric subtraction point the underlying one and two 
loop scalar master integrals involve structures not seen in the 3-point momentum configuration 
where one external momentum is nullified. In the latter case at low loop order one ordinarily 
only encounters rationals and the Riemann zeta function, C( z )- F° r the symmetric point one 
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finds the function 



s n {z) = i=3 




21 



for various arguments where Li n (z) is the polylogarithm function as well as one specific combi- 
nation of the harmonic polylogarithms which we denote by 



H 31 + H 43 



(2.10) 



The master integrals where these originally arise are summarized in [33] but the explicit evalua- 
tion are given in a set of articles, |44 [ I45 [ |4*6" 1 147| . Given that these functions will arise we record 
the relevant numerical values that we needed which are 



C(3) 
1 



1.20205690 , £ = 6.34517334 , tp' 



<■•'"(§) = 488.1838167 , s 2 [^J = 0.32225882 , s 2 
83 ' 2 



10.09559713 , 

= 0.22459602 



0.32948320 , s 3 - 
6 



= 0.19259341 

where ip(z) is the derivative of the logarithm of Euler T-function. 



(2.11) 



3 Renormalization. 



We devote this section to general aspects of MOM scheme renormalization. Having described 
the computer algebraic machinery used to construct the vertex functions we now recall how the 
MOM schemes are defined where we regard the amplitudes of (|2.4h as having been determined 
to the finite part. The divergences are removed into the coupling constant renormalization con- 
stant. However, to two loops this is an iterative procedure which is entwined with the 2-point 
function renormalization. This is because to extract the coupling constant counterterm from 
the vertex function one has to pay attention to the wave function renormalization constants of 
the external fields of the vertex function. Therefore, one first determines the one loop wave 
function renormalization constants in the MOM scheme of interest which is then fixed in ex- 
amining one loop vertex function defining that particular MOM scheme. For both the 2-point 
functions and the specific vertex the MOMi renormalization constant is defined so that at the 
subtraction point there are no 0(a) corrections after the renormalization constant is defined. 
Throughout we will use the syntax that in MOMi or equations i represents ggg, q or h. Once the 
one loop renormalization constants are fully determined then one repeats the exercise for the 
two loop contribution to first the 2-point function and then the associated vertex. The reason 
for explicitly defining the procedure is to ensure that there is no inconsistency in determining 
the coupling constant renormalization constant. The finite parts of the one loop wave function 
renormalization constants impact upon the finite parts of the two loop MOMi coupling con- 
stant renormalization constants, |21| . Otherwise an inconsistency in the renormalization group 
would emerge in trying to deduce the anomalous dimensions as well as the associated conver- 
sion functions for each renormalization constant. As a check on the vertex functions we have 
computed, we have verified that the two loop MS coupling constant renormalization constant 
of [H [21 [3l Hj correctly emerges when renormalizing at the symmetric subtraction point. One 
final point concerning our MOMi scheme renormalizations and that is that we first determine 
all the amplitudes before setting the renormalization constants for each scheme. For MOMi 
schemes we render that channel with the e divergences to have no 0(a) corrections. However, 
we emphasise that this is by no means the only way of defining the renormalization constants 
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within the MOM ethos. An alternative, for instance, is to first multiply the vertex function 
with its Lorentz tensor structure present, by the tensor of the corresponding Feynman rules. 
Then the coupling constant renormalization constant is defined by ensuring that there are no 
0(a) corrections to this object. We note that doing this for each of the three schemes would 
involve several of the non-Feynman rule amplitudes in contributing to the coupling constant 
renormalization. We have not proceeded in this way as it does not appear to be in keeping with 
[21|. I22j . However, in providing the full vertex structure in terms of the Lorentz tensors in the 
MS scheme an interested reader has the opportunity to study such alternative MOM scheme 
definitions. Indeed it may be the case that convergence of certain perturbative series could be 
improved in such a way. 

Having reviewed the procedure we followed we now comment on the relation of the param- 
eters of the theory in different schemes. For QCD the relevant parameters are the coupling 
constant and the linear covariant gauge fixing parameter. As the former is vertex dependent 
we comment on the latter first. As outlined the MOMi renormalizations require 2-point func- 
tion renormalization. Therefore, like the coupling constant, [21] , the gauge parameter can be 
different in different schemes. To relate them we follow the standard method and define 



«MOMi(^) 



7 MOMi 
7 MS 



o- 



MS 



(M) 



(3.1) 



where the subscript on the parameter refers to the scheme the variable is defined with respect 
to and Za is the gluon wave function renormalization constant. We follow the same conventions 
as [38] in defining the gauge parameter renormalization constant, Z a , by 



0'< 



Z a 1 Z A a 



(3.2) 



where the subscript, , indicates the bare parameter. With this convention then the respective 
anomalous dimensions satisfy, |48] . 



7a(a,a) 



7a(o,o) 



(3.3) 



in each scheme and j a (a, a) is the anomalous dimension of the linear covariant gauge parameter. 
In carrying out the renormalization in each of the three schemes we have determined OMOMi(/ i ) 
for each of the three cases and found that using the MOM scheme definition of [21] 



«MOMi 



1 + 

+ 



80T F N f - 9a|jg - 18a 



l MS 



97 



C A 



^MS 
36 



18^ - 18a|jg + 1900^ - 576C(3) aM g + 463^ + 864C(3) - 7143 
320a|^ + 3200^ - 2304C(3) - 4248] C A T F N f 



[4608C(3) - 5280] C F T F N f ] 



MS 
288 



+ 



O (c 



MS J 



a- 



MS 



(3.4) 



To ensure that the mapping is not divergent due to poles in e one has to iteratively solve (I3.ip 
order by order in perturbation theory hand in hand with the coupling constant of that scheme. 
This is because when we set the renormalization constants in a scheme both parameters of the 
explicit forms belong to that particular scheme. From (j3.4|) we see that the gauge parameter 
mapping is the same for all three schemes. Indeed it is the same as that for the RI' scheme, 
[48] . This is not unexpected as the parameter mapping is effectively the conversion function 
and reflects an underlying feature of the renormalization group. In essence it tracks how the 
scheme is defined for that one parameter amidst the renormalization of all the other parameters 
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and wave functions within the Green's functions. This will become evident later for other 
renormalizations. Therefore, given this feature the three loop term of (|3.4p has already been 
given in [H]. Further, we recall that it means that the Landau gauge is preserved between the 
schemes. 

The procedure to define the relation between the coupling constants in different schemes is 
similar. Though as we are dealing with three MOM schemes then the definitions are different 
in each case. We follow [21] and [22] for this and define 



a MOMses(/ i ) 



«MOMq(^) = a Ms(<") 



n MQMggg (p) 
u MOM q{p) ^ 



f OMq (p) N2 



J (l)MOMggg 



IlMS( p ) ( E MS(p) 



y-qqg / .,2 

(1)MS ( ~^ 



a MOMh(A i ) 



i-rMOMh 

9 



(p) (s 



MOMh 



(P) 



U MS (p) ( E MS(p) 



£ {l) S MOMq( 



(l)MS v 



J (l)MOMli( ^ ' ^ 



2 

(3.5) 



where we use the same designation for the vertices as before. Clearly the definitions involve 
the respective vertex functions evaluated at the symmetric subtraction point. Moreover, the 
amplitude chosen is that which has the divergences in e prior to renormalization or equivalently 
the amplitude which corresponds to the vertex Feynman rule. In the case of the triple gluon 
vertex we have chosen channel 1 which is only part of the Feynman rule. However, as will be 
apparent from the explicit results the other amplitudes from 2 to 6 are related in the way one 
would expect from the vertex structure so that our definition is consistent. Whilst in each of the 
three definitions the corresponding amplitude in the MOMi schemes have no 0(a) corrections, 
we have formally included it to ensure the normalization is correct and that the ratio of the 
vertex amplitudes from MOMi to MS begins with unity. The other main feature of (|3.5p is the 
presence of the 2-point functions. Specifically Ii g (p), T, c (p) and are respectively the scalar 

amplitudes of the gluon polarization and the Faddeev-Popov ghost and quark self-energies in the 
various schemes. The particular combination of which of these appears follows from the vertex 
of that MOMi scheme. In deriving the perturbative relations between these two parameters one 
has to proceed iteratively order by order in perturbation theory paying attention to the gauge 
parameter mapping in the same scheme at the same time. 

One particular property of the gauge parameter and coupling constant mapping between the 
MOMi and MS schemes is that we can now construct the other conversion functions for the wave 
function renormalizations and the coupling constant itself. Whilst the latter is not unrelated to 
(|3.5p we note that we regard conversion functions as being derived from the explicit forms of 
the renormalization constants themselves in the two schemes to be consistent with other work. 
Therefore, since we define the coupling constant renormalization constant, Z g , by, 
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then the conversion functions are given by 



a 



MOMi 



(a, a) 



Z 



MOMi 



Z MS 



a 



MOMi 



(a, a) 



Z 



MOMi 



z 



MS 



(3.6) 



(3.7) 



where 4> S {A, ip, c}. We will record the explicit forms of these for each of our schemes but note 
that to determine them we follow the same iterative procedure as we did in deriving (|3.4p . One 
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main benefit of the conversion functions and the parameter mappings is that we can deduce the 
/3-function and anomalous dimensions to three loops from the renormalization group. Specifically, 

122, 



P Ml ( a MOMi) «MOMi) 



and 



MOMi/ „ \ 

l<t> V a MOMi) a MOMiJ 



^ Q MOMi 



+ 



a Ms^ S ( a MS' a Ms) Q a 



<9 a MOMi 



MS 



MS->MOMi 
(3.8) 



^ ( a MsJ + P l a MsJ 







+ a 



MS 



MS / 
la \ 



a MS' a MS; f) n 



da 




lnC™.a, 



MS 



MS' a Ms) 



( a MS' a Ms) 



i r MOMi 

1U, MS' MS/ 



MS^MOMi 
(3.9) 



where the MOMi /3-functions will depend on the gauge parameter and only be scheme inde- 
pendent at one loop as these are mass dependent renormalization schemes, |21j . Though the 
MS /3-function is independent of a which is why it has only one argument. The mapping 
MS — > MOMi indicates that the object within the square brackets is first computed in terms 
of MS variables and then these variables are mapped back to the MOMi scheme variables by 
inverting (|3.4p and those derived from (|3.5p . As a check on this procedure we have calculated 
the /3-function and anomalous dimensions to two loops for each MOM scheme directly from the 
two loop renormalization constants and verified that the computation contained on the right 
hand sides of (|3.8p and (|3.9p are in agreement at the same order. This only requires the one loop 
terms of the conversion functions and therefore this provides a check on our computer routines 
designed to do this automatically. Finally, in order to ease comparison of our three loop MOMi 
results with the MS scheme in the same notation the numerical expressions for SU(3) for the 
latter scheme are, [H O O gj [5] , 



/3 MS (a) 



- [11.0000000 - 0.6666667 N f ]a 2 

- [102.0000000 - 12.6666667^/] a 3 

- 1428.5000000 - 279.6111111^/ + 6.0185185AT/ 2 
[0.6666667iV/ - 6.5000000 + 1.5000000a]a 



0(a 5 ) 



7 c MS (a,«) 



- 66.3750000 - 12.3750000a - 2.2500000a 2 - 10.1666667^/ 

- 915.9625108 - 165.2479023a - 33.9291631a 2 - 5.9062500a 3 

- [186.8599000 - 9.0000000a] N f + 7. 9629630 A^/j a 3 + 0(a 4 

[0.7500000a - 2.2500000]a 

- [17.8125000 + 0.5625000a - 1.2500000TV/] a 2 

- 256.2687446 - 2.1729239a + 0.5114565a 2 - 1.2656250a 3 

- [46.2756056 - 3.9375000a] N f - 0.9722222A/ 2 ] a 3 + 0(a 4 ) 
1.3333333aa 



+ 



22.3333333 + 8.0000000a + l.OOOOOOOa 2 



1.3333333 A, 



528.3243079 + 109.4435121a + 20.0342561a z + 3.7500000a J 

- [61.1111111 + 8.5000000a] A/ + 0. 74074077V/ 2 1 a 3 + 0(a 4 ) (3.10) 
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to the same numerical accuracy. To ease the presentation on the eye we will use the convention 
that when the scheme appears on the function on the left hand side then the variables, such as 
the coupling constant and gauge parameter, are the variables in the same scheme. The exception 
to this is the conversion functions, Cf om (a,a) and Cf OMi (a,a), where the arguments are the 
MS variables. 



4 Triple gluon vertex. 



We now begin the mundane task of recording our results for each scheme and its respective 
vertices by concentrating on the triple gluon vertex first. As indicated earlier the full analytic 
versions of all the results for an arbitrary gauge, in this and the next two sections, have been 
included in a separate data file. The MS SU (3) numerical values for the amplitudes are 



MS 



MS 



J (8) 



MS 



MS 



J (2) 

2 ( 5 ) 
- 1 - 



MS 

(p,q) 



2 "( 3 ) 

= - ^ 



MS -(6)(P>9) 
1.1212444 - 3.7618956a - 
- 0.04173667V)] a 



MS 



MS 



MS 



1.2890232a 2 + 0.1250000a 3 



+ 



29.7530676 + 16.4600770a - 9.7794300a 2 - 3.2060809a 3 
- 1.6522848a 4 + 0.2812500a 5 



- [11.5677203 - 0.9686976a - 0.9112399a 2 + 0.4166667a 3 ]iV / 
+ 0(a 3 ) 



MS " ^ 14 ) 



Mms 



+ 



0567163 - 3.3280464a - 0.5079304a 2 + 0.0573179a 3 - 1.0926858iV)] a 
116.0789643 - 13.6830818a + 0.3484134a 2 + 4.7763124a 3 



+ 0.8908609a 4 - 0.1289652a 5 



- [20.2710109 + 1.0153018a - 0.5745217a 2 - O.igiOSgea 3 ]^/ 
+ 0(a 3 ) 

7.3683002 - 3.3518377a - 0.5701159a 2 + 0.1926821a 3 - 1.21300967V) 



+ 



126.0048710 - 11.8048854a + 3.7795690a 2 + 4.3779190a 3 
+ 1.2887087a 4 - 0.4335348a 5 



- [23.5898191 - 0.0155813a - 0.9363317a 2 - 0.6422738a 3 ]iV / 
+ 0(a 3 ) 

- 0.3115839 - 0.0237913a - 0.0621855a 2 + 0.1353643a 3 - 0.12032387^ 

- 9.9259067 + 1.8781964a + 3.4311557a 2 - 0.3983934a 3 
+ 0.3978478a 4 - 0.3045696a 5 

- [3.3188082 - 1.0308831a - 0.3618100a 2 - 0.4512142a 3 ]7V/ 
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+ 0(a 3 ) 



(4.11) 



We have indicated the relations between amplitudes of the various projection tensor channels. 
These are consistent with the expectations for the structure of the vertex from symmetry given 
that we have evaluated the vertex function at the symmetric point. In addition in this context 



\7) 



(p, q) 



MS 



J (8) 



\p, q) 



MS 



J (10) 



(P. 9) 



MS 



(4.12) 



which we checked was true to two loops prior to the numerical evaluation. In Table 1 we give 
the comparison of our results for the channel 1 amplitudes to those of [22] . In this and other 
comparisons we have used the coupling constant convention of the presentation in [22] where 
the series was in powers of a s = g 2 /(A7r 2 ). This is to retain the error estimates given in [22] , It 
is clear that our evaluation is comfortably close to the approximations of 



s (l s )>></) 


MS 


r 2 


C A T F N f 


C F T F N f 


Ref [22] 
This paper 


0.22(4) 
0.2066185 


- 0.65(7) 
- 0.6620808 


0.408(10) 
0.4052081 



Table 1. Comparison of channel 1 two loop Landau gauge amplitude with [22J by colour factor. 



The MOMggg scheme amplitudes satisfy the same relations and in particular we have 

1 

MOMggg 
J (4)'(^9) 



MOMggg 



J (7) 



'(p,q) 



MOMe 



J (8) 



'(p,q) 



MOMe 



J (2) 

- s 

- s 



(p,q) 



vest 
2 ( 3 ) 
1 



Xp,q) 



(6) 



'(j>,q) 



2£§f(p )9 ) 



MOMggg 
MOMggg 
MOMrkk 



2 ( 5 ) 



MOMggg 
MOMggg 



- 1 + 0(a 3 ) 



MOMe 



J (14) 



(p,q) 



MOMggg 



7.0567163 - 3.3280464a - 0.5079304a 2 + 0.0573179a 3 
- 1.09268587V/] a 

- 78.7833165 - 99.1996625a + 10.0012225a 2 + 10.9109237a 3 

- 1.2024953a 4 - 0.2831609a 5 + 0.0214942a 6 

- [34.3080792 - 16.2422875a - 1.8203921a 2 + 0. 6079935a 3 ]N f 
+ 3.7791007^/1 a 2 + 0(a 3 ) 



J (13) 



(p,q) 



MOMggg 



7.3683002 - 3.3518377a - 0.5701159a 2 + 0.1926821a 3 
- 1.2130096^/] a 

- 77.4614044 - 103.6568232a + 5.5514994a 2 + 12.5463340a 3 

- 2.9431069a 4 - 0.5253739a 5 + 0.07225580a 6 

- [35.3894861 - 16.0837324a - 1.7300353a 2 + 1.1212780a 3 ]iV / 



+ 4.1952457A/ 



a 2 + 0(a 3 
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Eg?)(p,g) 



MOMggg 



S flf)(P^) MQMggg 

0.3115839 - 0.0237913a - 0.0621855a 2 + 0.1353643a 3 
- 0.1203238^/] a 

1.3219120 + 4.4571607a + 4.44972304a 2 - 1.6354103a 3 

+ 1.7406116a 4 + 0.2422130a 5 - 0.0507616a 6 

+ [1.0814069 + 0.1585516a + 0.0903568a 2 + 0.5132845a 3 ]iV) 

a 2 + 0(a 3 ) (4.13) 



- 0.4161497A/ 



with the corresponding relation 
^(7?^) MOMggg 



J (8) 



MOMggg 



J (io) 



(P,9) 



MOMggg 



(4.14) 



also being true to two loops analytically. Given the nature of the MOMggg scheme the relations 
for the amplitudes of channels 1 to 6 demonstrate that our renormalization is consistent and 
that our projection has been implemented consistently within our Form programmes as well as 
being a check on our Reduze database. Next we record the four conversion functions for the 
wave function and coupling constant renormalizations in the MOMggg scheme relative to the 
MS scheme. These are 



c MOMggg (0)a) 



13.2462444 - 1.5118956a - 0.1640232a 2 



+ 0.1250000a 3 - 1.7084032^ 
217.0368707 + 36.7247782a + 7.0535877a 2 - 1.1557619a 3 
+ 1.0453915a 4 - 0.0996192a 5 - 0.0156250a 6 



^ OMggg (a,a) = 



c M0Mggg (OjQ) 



C^ 0Mggg (a,a) 



- [33.1527255 + 3.7086335a - 0.0047429a 2 - 0.6354341a 3 ]iV / 

- 0.5342658A/1 a 2 + 0(a 3 ) 

1 + [8.0833333 + 1.5000000a + 0.7500000a 2 - l.llllllliV/1 a 
+ [256.1034914 + 31.4182743a + 8.4375000a 2 + 2.8125000a 3 

- [53.9129277 + 1. 6666667a] iV> + 1. 23456797V/ 1 a 2 + 0(a 3 ) 
1 + 3.0000000a 

+ [80.9357699 + 3.0725670a + 1.3465290a 2 - 5.9375000^/ 
+ 0(a 3 ) 

1 - 1.3333333aa 

25.4642061 + 11.5753171V, 2. 7222222a 
+ 0(a 3 ) 

MOMe 



1 2.3333333iY/J a 2 



(4.15) 



where we note that the expressions for all bar Cf UMggg (a , a) are formally the same as those 
for the RI' scheme. By definition C^ 1 (a, a) is unity. From C g ggg (a, a) we can deduce the 
explicit form of the coupling constant mapping. Though as the full form is cumbersome we 
record the analytic version for the Landau gauge which is 



«MOMggg 



a MS + 



-46tt 2 + 1188 



C A + 



128vr z - 192y/(±) - 432 T F N f 



MS 
162 



13 



+ 



19044ft// (|)) 2 - 25392ttV(§) - 6938784<//(|) - 100602^ w (i) 

- 72643392s 2 (§) + 145286784s 2 (§) + 121072320s 3 (§) - 96857856s 3 (§ ; 
+ 2767367T 4 + 46258567T 2 - 113724S + 8301852C(3) + 40126833 

- 504468^^ + 6053616^ + 541836^ 

V3 v3 v3 



+ 



+ 



+ 



141312ttV(±) - 105984(V' / (i)) 2 - 2960064^//(±) + 33592320s 2 (f ) 
- 67184640s 2 (§) - 55987200s 3 (f ) + 44789760s 3 (§) - 47104vr 4 



19733767T 2 + 2239488S - 8957952C(3) 



ln 2 (3)vr ln(3)vr vr 
+ 233280 y -J 2799360^= 250560^ 



26695008 

.3 



C A T F N f 



\/3 y/3 
124416^"(|) - 1492992^' (§) - 331776tt 4 + 995328tt 2 
- 4478976S + 6718464C(3) - 7138368] C F T F N f 
147456(V>'(|)) 2 - 196608ttV(|) + 2322432^(1) + 65536tt 4 



1548288vr 2 + 2923776] T F Nj 



7 3 
MS 



419904 



+ 



°( a m) 



By contrast the full gauge dependent version in numerical form for SU (3) is 



(4.16) 



OMOMggg 



*MS 



+ 



+ 



26.4924889 - 3.02379130^ - 0.32804640^ 



+ 0.25000000^ - 3.41680647V} 



2 

TZS 



960.4627167 - 46.71207940^ + 7.92851320^ + 9.1110752o| 



MS 



MS 



+ 1.03757210^ 



0.3222558o|jg + 0.0156250o|^ 



202.0850109 - 8.0802973a 



MS 



1.6907923a 



MS 



+ 0.010434101^] Nf + 7.6873930^/] a 3 ^ + O (a 4 ^) . (4.17) 



Equipped with these we have computed the three loop renormalization group functions for the 
MOMggg scheme. Again given space considerations we record the Landau gauge expressions for 
each case. They are 



MOMggg 



(o,0) = - 



+ 



l ^C A -\T F N f 



a 2 - 



y C\ - AC F T F N f - ^C A T F N } 



[ [209484^' {\)f - 279312vrV( 1 r) + 37087200^(1) 
+ 368874t//"( 1 t) + 266359104s 2 (£) - 532718208s 2 (§) 

- 443931840s 3 (f ) + 355145472s 3 (f ) - 890560vr 4 

- 247248007T 2 + 416988S - 30440124C(3) - 51650217 



+ 1849716 



ln 2 (3)7T 



22196592 



ln(3)7T 

~7T 



1986732 



7T 



c\ 



+ 



1656000vr 2 V' / (5) 



1242000(V>' ( \ )f 



38988864V>\ 3 ; 
134136V/"(§) - 220029696s 2 (|) + 440059392s 2 (f ) 

4 



+ 366716160s 3 (f ) - 293372928s 3 (f ) - 194304vr 4 
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+ 25992576vr 2 - 8363088S + 43914960C(3) + 49845132 



1527984 



ln 2 (3)vr 



+ 18335808 



ln(3)vr 



7T 



+ 1641168-= 

V3 



+ 



3- 



y/3 \/3 
2045952(/(i)) 2 - 27279367rV( i f ) + 11591424-0' ( 
+ 44789760s 2 (f ) - 89579520s 2 (§) - 74649600s 3 (f 
+ 59719680s 3 (§) + 909312^ 4 - 7727616tt 2 
+ 2985984S - 11943936C(3) - 8460288 



C\T F N f 



+ 311040 



ln 2 (3)7r 



- 3732480 



ln(3)7r 



7T 



- 334080—= 
V3 



C A T F N? 



+ 



+ 



+ 



>/3 ^3 
786432ttV(§) - 589824(V>'( ±)) 2 - 442368^'(|) 

- 262144vr 4 + 294912vr 2 - 82944] TpNf 

4758912^'(i) - 456192^'"(i) + 1216512tt 4 - 3172608tt 2 
+ 16422912S - 24634368C(3) + 23421312] C A C F T F N f 
165888^"' (£) - 4423687T 4 - 5971968S + 8957952C(3) 

- 7091712] C F T F Nf - 839808C| TpiV/] 



419904 



+ 0(a e 



(4.18) 



- [13CU - 8T F N f ] - 



+ 



C\ + 7776C F T F N f 



+ 



+ 



+ 



1794y/(±) - 1196vr 2 - 2655 
4064vr 2 - 6096^(1) + 2304 C A T F N f 

3072V>'(|) - 2048vr 2 + 1152 T F Nj 

2310672^V(±) - 1733004(^'(i)) 2 
- 1307826^"' (±) - 944364096s 2 (f ) + 1888728192s 2 (§) 
+ 1573940160s 3 (f ) - 1259152128s 3 (§) + 2717312vr 4 
+ 80915616tt 2 - 1478412S + 164768580C(3) - 117299583 



1944 

121373424?// (§) 



ln 2 (3 vr In 3)vr vr 3 

6558084 \J- + 78697008 V + 7043868—= 

V3 v/3 V3 



+ [l071108(V/(±)) 2 - 14281344vrV( 1 r) + 134602560^(1) 
+ 804816<'(i) + 1017847296s 2 (f ) - 2035694592s 2 (f ) 

- 1696412160s 3 (f ) + 1357129728s 3 (|) + 2614272^ 4 

- 89735040^ 2 + 30023136S - 100567008C(3) + 80188056 



+ 7068384 



l n 2 (3)7T 

~7T 



- 84820608 



ln(3)vr 

7T 



7T° 



- 7591968^= 
\/3 



C\T F N } 



+ 



25804800vrV (§) - 19353600^' (|)) 2 - 40849920^(1) 

- 268738560s 2 (f) + 537477120s 2 (|) 

+ 447897600s 3 (f ) - 358318080s 3 (§) - 8601600vr 4 

+ 27233280vr 2 - 17915904S + 17915904C(3) + 2052864 



15 



- 1866240 



ln 2 (3)7T 

~7T 



+ 22394880 



^fc + 2004480^ 



C A T 2 F Nf 



+ 



+ 



+ 



8257536(?//(±)) 2 - 11010048ttV (§) + 6193152^ 



'<\) + 36700167T 4 



4128768tt 2 + 1161216 



T F Nf 



1617408<'(§) - 27993600^(|) - 4313088tt 4 + 18662400tt 2 
- 58226688S - 147806208C(3) + 125971200] C A C F T F N f 
35831808?// (i) - 995328<'(±) + 2654208tt 4 - 23887872vr 2 
+ 35831808S + 53747712C(3) - 47029248] C F T F Nf 



- 5038848C F T F N f 



a 



2519424 



+ 0(a 4 ) 



(4.19) 



MOMggg 



(a,0) 



-C A a 



+ 



+ 



138^'d) - 92vr 2 - 63 



+ 



C A 



256vr 2 - 384/(i) + 72 T F N f 



1 C A a 2 
432 



177744vrV(±) - 133308^' (^)) 2 - 9492336^' (£) 
- 100602?//" (§) - 72643392s 2 (f ) + 145286784s 2 (§) 
+ 121072320s 3 (f ) - 96857856s 3 (§) + 209024vr 4 
+ 6328224tt 2 - 113724S + 11993508C(3) - 9641835 



504468 



+ 



ln 2 (3)7r 



+ 6053616 



ln(3)vr vr 3 
— M~ + 541836^ 

V3 y/3 



741888(?//(±)) 2 - 989184vrV(±) + 5019840?//(±) 
+ 33592320s 2 (f ) - 67184640s 2 (§) 

- 55987200s 3 (f ) + 44789760s 3 (f ) + 329728vr 4 

- 3346560vr 2 + 2239488S - 5318784C(3) + 6350400 

3 



ln 2 (3)vr ln(3)vr 
+ 233280 W 2799360 



250560-^= 



+ 



+ 



1376256vrV(i; 



1032192(V' / (|)) 2 



- 458752vr 4 + 73728vr 2 



124416V>'(|) - 1492992?// ( 



V3 

110592V? 

1762560 

331776vr 4 



C A T F N f 



1 (V 

3' 



T F Nf 



+ 995328^ 2 - 4478976S + 1399680 



C F T F Nf 



C A a 3 
559872 



+ 0(a 4 



(4.20) 



and 



MOMggg, n s 
7^ sss (a,0) 



[25C A - 6C F - 8T F N f ] 



C F a 2 



+ 



4600tt 2 - 6900?/> 



3' 



39690C(3) + 61011 



r 2 



+ 



1656^'d) - 1104vr 2 + 15552C(3) - 23760 



CaCf 



16 



+ 
+ 



21408^/(i) - 14272tt 2 + 10368C(3) - 28800 C A T F N f 



+ 

+ 0(a 4 



3072vr 2 -4608^(1) -4320 C F T F N f + 1944C| 

C F a 2 



4096vr 2 - 6144^'d) + 1152 



T F Nf 



1296 



(4.21) 



That for the quark anomalous dimension is relatively compact as there is no one loop term in the 
Landau gauge. The explicit gauge dependence is given in the numerical evaluations for SU(3) 
which are 



MOMggg 



(a, a) 



MOMg 

7a 



'{a, a) 



MOMggg 



(a, a) 



[11.0000000 - 0.6666667iV)]a 2 
102.0000000 + 19.6546434a - 0.2710840a 2 - 5.8591391a 3 
+ 1.1250000a 4 

- 12.6666667 + 2.0158609a + 0.4373952a 2 - 0.5000000a 3 ] N f 
1570.9843804 + 658.0709292a + 269.2238338a 2 + 43.0029610a 3 

- 99.2797189a 4 + 14.8550247a 5 + 5.3345924a 6 - 0.7031250a 7 



+ 



0.5659290 - 43.2393672a - 22.7471960a 2 - 19.8709555a 3 



N f 



+ 14.8347569a 4 + 0.9764184a 5 - 0.2812500a 6 

67.0895364 + 4.6479610a + 0.8898051a 2 - 2.3056953a 3 
3l „4 , r>,^ 



Nf 



+ 2.65811557V/ a 4 + 0(a 5 ) 
[0.6666667iY/ - 6.5000000 + 1.5000000a]a 



+ 



16.9095110 - 41.6433767a + 6.1533855a 2 + 0.9920696a 3 

- 0.3750000a 4 

12.0931233 - 5.4744039a + 0.2813024a 2 + 0.1666667a 3 



+ 1.53713027V/ 

1308.9386744 - 647.9260677a + 376.2301295a 2 + 6.3971133a 3 



33.0162468a 4 + 7.3253130a 5 + 1.0008734a 6 



0.1640625a 



7 



491.4309500 - 302.3530495a + 52.3029150a 2 + 6.3604344a 3 



N 



/ 



- 6.7153148a 4 - 0.1288604a 5 + 0.0729167a 6 
74.9190172 - 29.3999306a + 1.4158902a 2 + 1.3449889a 3 
- 6. 20226947V/ ] a 3 + 0(a 4 ) 
[0.7500000a - 2.2500000]a 



+ 



N f 



+ 



8.7955100 - 21.1728971a + 2.6547391a 2 + 1.3710348a 3 
- 0.1875000a 4 - [4.4378145 - 1.7292715a] TV/ 
548.8492387 - 436.6720556a + 199.2938032a 2 + 32.7086138a 3 



- 30.0123943a 4 + 1.4303427a 5 + 0.9535617a 6 



0.0820312a' 



157.4669179 - 127.5457558a + 20.0522192a 2 + 7.7592757a 3 
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MOMggg, N 



- 1.5131126a 4 
1.3333333aa 



Nf + [19.9741163 - 6.9775531a] Nf a 3 + 0(a 4 ) 



+ 



22.3333333 - 10.5455407a + 9.0317217a 2 + 1.4373952a 3 
- 0.3333333a 4 - [1.3333333 — 3.0742604a] Nf 
94.7943290 - 204.1998798a + 218.8404110a 2 - 30.4216658a 3 



34.4073860a 4 + 6.3159939a 5 + 1.2577208a 6 



0. 1458333a ' 



76.8672720 - 80.5601965a + 53.0718118a 2 + 7.0576676a 3 



2.6899779a 4 



N f + [5.2596320 - 12.4045389a] Nf 



+ 0(a 4 ) . 



(4.22) 



We recall, |22| . that unlike the MS scheme the three loop term of the MOMggg /3-function is 
cubic in Nf and not quadratic. 



5 Ghost-gluon vertex. 



We repeat this exercise now for the structure of the ghost-gluon vertex and the associated MOMh 
renormalization scheme. Though there are fewer amplitudes than for the triple gluon vertex. 
For the MS scheme the two amplitudes are 



MS ' 



2.2324710 + 0.3280464a - 0.1464920a 2 



MS 



49.2999213 + 16.6398011a + 2.3538026a 2 + 0.0885284a 3 
+ 0.1098707a 4 - [4.0701546 + 0.4453256a + 0.1627713a 2 ]iV / 
+ 0(a 3 ) 

1.4824710 - 0.1640232a - 0.1464920a 2 
+ 



35.1253580 + 2.2852097a - 0.3277110a 2 + 0.0885284a 3 
+ 0.1098707a 4 - [3.2368212 + 0.4453256a + 0.1627713a 2 ]iV / 
+ 0(a 3 ) 

leading to the MOMh scheme expressions 



a 

(5.23) 



J (2) 



MOMh 
MOMh 



- 1 + 0(a 3 ) 

1.4824710 - 0.1640232a - 0.1464920a 2 



+ 



4.3185039 + 0.6852156a + 0.0493136a 2 - 0.0591276a 3 



- 0.0643817a 4 
+ 0(a 3 ) . 



[1.5896312 + 0.4453256a]A^ 



(5.24) 



As channel 1 contained the poles in e after MOMh renormalization then in that scheme there 
are no corrections at the symmetric subtraction point. Unlike the other two vertices we do 
not include a table comparing our numerical expressions with those of [22] as there is no direct 
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comparisoiill. Moreover, we have examined whether various combinations of our amplitudes 
can produce agreement, because of potentially different conventions of defining the basis, but 
have failed to find any. For other quantities related to the MOMh scheme such as the coupling 
constant mapping and some coefficients of the three loop SI/ (3) Landau gauge /3- function we 
find similar but minor discrepancies with [22] which suggest a consistent typographical error in 
the presentation of certain equations in [22]. Further, related inconsistencies for the results for 
this vertex given in [22] have also been noted in |28j. Whilst we will comment in more depth 
in section 7 in the context of the results of the other two schemes we note that our benchmark 
check, |24j . on the Landau gauge SU(3) f3- function for MOMh was slighty better than that of 
the MOMggg case. 

Continuing with the presentation of our results, the SU (3) numerical values for the respective 
MOMh conversion functions to MS are 
-(MOMh/ 



C, 



l (a, a) 



1 



9.2741377 + 1.0780464a + 0.2285058a 2 - 0.5555556^/ 
191.9555891 + 18.6287007a + 1.8982822a 2 + 0.7339954a 3 
- 0.0675921a 4 - [27.3210789 - 0.1535891a - 0. 3808430a 2 ]N f 



0.1543210^ 



+ 0(a 3 ) 



^yMOMh 



(a, a) 



1 + 



8.0833333 + 1.5000000a + 0.7500000a 2 - 1.111111L/Vj 



+ 



256.1034914 + 31.4182743a + 8.4375000a 2 + 2.8125000a 3 



C, 



MOMh 



(a, a) 



- [53.9129277 - 1.6666667a] A) + 1.23456797V/ 
1 + 3.0000000a 



+ 0(a 3 ) 



+ 



5.9375000/V/ 



^MOMh 

V 



(a, a) 



80.9357699 + 3.0725670a + 1.3465290a 
+ 0(a 3 ) 

1 — 1.3333333aa 
- 25.4642061 + 11.5753172a + 2.7222222a 2 - 2.33333337V) 
+ 0(a 3 ) . 

The analytic form of the first of these leads to the coupling constant mapping 



(5.25) 



°MOMh 



+ 



15-0'd) - 10vr 2 + 6151 C A 



MS 
108 



+ 



24QT F N f 

450(^'(i)) 2 - 600vrV(§) - 458928^' (§) - 3213^'"( 1 r) 
- 3825792s 2 (§) + 7651584s 2 (f ) + 6376320s 3 (f ) - 5101056s 3 (f ; 
+ 8768vr 4 + 305952vr 2 + 7776S + 153576C(3) + 6521760 



- 26568^^ + 318816^ + 28536^ 



C\ + 460800TIAT/ 



+ [206784^(|) + 1492992s 2 (f ) - 2985984s 2 (f ) - 2488320s 3 (|) 
+ 1990656s 3 (f ) - 137856vr 2 - 995328C(3) - 4015296 



+ 10368^2^ - 124416- ln(3) " 



V3 



V3 



11136 



7T 



V3 



C A T F N f 



[1492992C(3) - 1710720] C F T F N f ] 



a- 



MS 



93312 



+ 



O (a 



MS) 



(5.26) 



'We are grateful to Prof. Chetyrkin for comments on this point. 
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or numerically, for an arbitrary linear covariant gauge, 
a MOMh 



18.5482754 + 2.15609280^ + 0.45701160^ 



l.llllllliVf 



+ 



641.9400674 + 97.24510470^ + 19.99828180^ + 2.94602990^ 



+ 0.02146060^ 



+ 1.2345679 AT/ 



a MS + 



85.5559502 + 3.28630980^ N f 
(4s) • 



(5.27) 



Using the renormalization group the three loop renormalization group functions emerge. Ana- 
lytically in the Landau gauge we have 



/3 MOMh (a,0) = - 
+ 



a 2 - 



y C\ - 4C F T F N f - ™C A T F N f 



4950<y (i)) 2 - 6600vr^'(i) + 2370816^'(i) 

+ 11781V>'"(|) + 14027904s 2 (|) - 28055808s 2 (§) - 23379840s 3 (f ) 
+ 18703872s 3 (f ) - 29216tt 4 - 1580544vr 2 - 28512S - 563112C(3) 

r 3 



11733336 + 97416±^ - 1168992 l ^ 



104632^= 



+ 



V3 

2400vr 2 ^'(i) - 1800(^»'(i)) 2 - 1864512^(1) - 4284V> W ( ±) 

- 10575360s 2 (f ) + 21150720s 2 (§) + 17625600s 3 (f ) 

- 14100480s 3 (f ) + 10624vr 4 + 1243008vr 2 + 10368S + 3854304C(3) 



+ 9722592 - 73440^4=^ + 881280^2^ + 78880-^jL 
V3 \/3 VS 



C 2 A T F N f 



+ [352512V/ (±) + 1990656s 2 (f ) - 3981312s 2 (f ) - 3317760s 3 (f ) 
+ 2654208s 3 (f ) - 235008vr 2 - 1327104C(3) - 1368576 



+ 13824 



ln 2 (3)vr In 3)vr vr 3 

\d- - 165888^^ - 14848^ 

V3 



C A T 2 F Nf 



+ 



34560tt 2 - 51840^' (i) - 5474304C(3) + 627264o] C A C F T F N f 

4 



+ [1990656C(3) - 1907712] C F T F Nf - 186624C F T F N f 
+ 0(a 5 ) 



a 



93312 



(5.28) 



. MOMh 

lA 



(a,0) = - [13Ca — 8TpNj] — 



+ 



+ 



195^(1) - 130vr 2 - 3186] C\ + 2592C F T F N f 



+ 



807T' 



120^(1) + 2808 



C A T F N f 



648 



54600vrV(§) - 40950(V>'(±)) 2 - 7120224^'(|) - 41769V>'"(±) 

- 49735296s 2 (f ) + 99470592s 2 (f ) + 82892160s 3 (f ) - 66313728s 3 (§) 

+ 93184tt 4 + 47468167T 2 + 101088S + 14628600C(3) - 37070136 



345384 



ln 2 (3)vr 



+ 1414608 



ln(3)vr 



7T*- 



+ 370968^= 
V3 



+ 



V3 V3 
25200(^'(|)) 2 - 33600vrV(|) + 7615296^(1) + 25704^ 



in i i ~ 

3- 
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+ 50015232s 2 (f ) - 100030464s 2 (§ ) - 83358720s 3 (f ) 

,ln 2 (3)7r 



+ 66686976s 3 (§) - 57344vr 4 - 5076864vr 2 - 62208E 



+ 4121280C(3) + 35848360 + 347328- 



- 4167936 



ln(3)vr vr 3 
— M~ - 373056^= 
V3 V3 



C\T F N f - 1119744C|T F iV / 



+ [23887872s 2 (f ) - 2115072^' (±) - 11943936s 2 (f ) + 19906560s 3 (§) 
- 15925248s 3 (§) + 1410048tt 2 - 3981312C(3) - 5225472 



- 82944 



ln 2 (3)vr 



+ 995328^^+89088^ 



C A T 2 F Nf 



+ 



414720vr 2 - 622080^' (|) - 32845824C(3) + 33716736 C A C F T F N f 



+ [11943936C(3) - 11446272] C F T F N? 



559872 



+ 0(a 4 ) 



(5.29) 



MOMh 

7c 



(a,0) 



and 



.MOMh 

lip 



(a,0) 



-C A a + \\lb<4>'(\) - 10vr 2 - 198 



C A + 72T F N f 



C A a 2 
144 



+ 



4200vrV(i) - 3150(^'(i)) 2 - 559008V>'(±) - 3213^"'(§) 

- 3825792s 2 (f) + 7651584s 2 (f ) + 6376320s 3 (f ) - 5101056s 3 (§) 

+ 7168vr 4 + 3726727T 2 + 7776S + 973944C(3) - 2855736 



26568 



ln 2 (3)7r 



+ 318816 



>"+t + 28536 +l 



+ [247104^'(i) + 1492992s 2 (f ) - 2985984s 2 (f ) - 2488320s 3 (f ) 
+ 1990656s 3 (f ) - 164736vr 2 - 186624C(3) + 2260224 

r3 



+ 10368^ - 124416- ln(3)7r 



V3 



- 4147207^/ + 186624C F T F N f 



- 11136^= 

V3 

C A a 3 



C A T F N f 



124416 



+ 0(a 4 



(5.30) 



[25C A - 6C F - 8T F N f ] 



C F a 2 



+ 



500vr 2 - 750V>'(|) - 13230C(3) + 29187 



+ 
+ 



180y/(|) - 120vr 2 + 5184C(3) - 10044 C A C F 
240^(|) - 160^ 2 + 3456C(3) - 14832] C A T F N f 



- 864C F T F N f + 648Cf + \lh2T F Nf 



C F a 2 
432 



+ 0(a 4 ) . (5.31) 



For comparison with MS and the other two MOM schemes we note the full arbitrary linear 
gauge values are 



P 



MOMh 



(a, a) 



11.0000000 - 0.6666667^/ }a 2 

102.0000000 - 14.0146029a - 2.7070116a 2 + 1.3710348a 3 



12.6666667 - 1.4373952a - 0.6093488a 2 



N f 
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MOMh 

lA 



(a, a) 



MOMh 

IC 



(a, a) 



2813.4929483 - 138.6933593a - 35.8596355a 2 + 4.4867569a 3 

+ 9.6466346a 4 + 1.3338552a 5 

617.6471542 - 40.9610609a - 14.8518881a 2 - 2.2357888a 3 



- 0.2607456a 4 N f + 21.5028181^/ a 4 + 0(a 5 ) 

[0.6666667^/ - 6.5000000 + 1.5000000a]a 

- 34.7278768 - 3.9421899a - 3.4864362a 2 + 1.8105174a 3 

8.1900386 - 1.4373952a - 0.8046744a 2 ] N f 
1195.1013833 - 28.3678253a - 8.6742911a 2 - 5.5422477a 3 
+ 10.3727415a 4 + 1.0931120a 5 

- 323.9161114 - 37.7892337a - 11.8481326a 2 - 1.7973717a 3 

- 0.1108276a 4 ] N f + [12.2638016 - 0.5937675a] iV/] a 3 + 0(a 4 
[0.7500000a - 2.2500000]a 

- 9.0788804 + 3.5599978a - 0.5362065a 2 - 0.2197413a 3 

- 0.7500000iV)] a 2 

- 462.6953814 - 89.8527064a - 1.2245066a 2 + 2.9101869a 3 

- 0.3430542a 4 - 0.1126679a 5 

- 76.4163746 - 1.4615295a - 0.0175290a 2 

+ 0(a 4 ) 
1.3333333aa 

22.3333333 + 0.0467440a + 2.1252097a 2 + 0.3906512a 3 

- 1.33333337V/] a 2 

+ r 



N f + 2.5000000^/ 



+ 



260.0472082 - 162.9606897a - 38.3957984a 2 + 2.9734643a 3 

Nf + 0.8888889 Nf 



+ 1.2107696a 4 + 0.2002985a 5 



- 47.3050219 - 12.9383858a - 2.4062325a 2 
+ 0(a 4 ) . 

The /3-function has the same Nf polynomial structure as the MS scheme. 



a 

(5.32) 



6 Quark-gluon vertex. 

Finally, we complete our task with the quark-gluon vertex and the MOMq scheme expressions. 
The MS amplitudes are 



MS 



= 1 + 4.3162206 - 0.5887601a - 0.4570116a 2 

+ [89.2876778 - 2.5488660a + 0.7959457a 2 + 0.2344278a 3 

+ 0.3427587a 4 - [12.1366772 + 0.9766280a + 0.5077907a 2 ]iV) 
+ 0(a 3 ) 
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£$) V (p, q) 



MS 



MS 



MS 



MS 



2.5980335 - 2.3056953a - 0.4140232a 2 

26.4812470 - 21.7488508a - 5.3984938a 2 + 0.4547874a 3 
+ 0.3105174a 4 - [6.2718940 + 1.0339459a + 0.4600258a 2 ]iV / 
0(a 3 ) 



MS 



2.0502686 - 2.5226306a - 0.5000000a 2 



+ 



12.7352941 - 25.2299763a - 6.6819786a 2 + 0.0320681a 3 
+ 0.3750000a 4 - [4.8715927 + 0.9193101a + 0.5555555a 2 ]iV / 
+ 0(a 3 ) 

4.3622718 + 2.3439072a + 0.5859768a 2 ] a 
131.9911153 + 45.4675027a + 4.8573516a 2 + 1.2207850a 3 



0(a 



0.4394826a y 

3\ 



[10.9228497 + 1.9532560a - 0.65108530^/ 



a 

(6.33) 



where the symmetry of the exchange of the two external quark legs is manifest. This emerged 
from the computation and was not imposed. In Table 2 we have given the comparison table 
analogous to Table 1 for this vertex. Again the numerical estimates and the results we have 
produced are virtually identical. 





MS 


r 2 


CfCa 


r 2 


C A T F N f 


C F T F N f 


Ref [22] 
This paper 


0.206(4) 
0.2048069 


- 0.20(4) 
- 0.2158263 


0.679(1) 
0.6755204 


- 0.4968(4) 

- 0.4958508 


- 0.0211(4) 

- 0.0221492 



Table 2. Comparison of channel 1 two loop Landau gauge amplitude with [22] by colour factor 



The corresponding MOMq scheme expressions are 



J (2) 



MOMq 
MOMq 



MOMq 



1 + 0(a 3 



MOMq 



2.5980335 - 2.3056953a - 0.4140232a z 



28.1605810 - 15.7267125a + 11.9916906a 2 + 5.1627786a 3 



+ 0.5676402a 4 + [3.3851901 + 1. 0339459a) N f 



a 2 + 0(a 3 



MOMq 



2.0502686 - 2.5226306a - 0.5000000a 2 

30.3859453 - 13.4480431a + 14.1587135a 2 + 6.3930196a 3 
+ 0.6855174a 4 + [2.5935165 + 0.9193101a)iV/l a 2 + 0(a 3 ) 
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MOMq 



4.3622718 + 2.3439072a + 0.5859768a 2 
40.2438356 + 27.9113517a + 15.1059206a 2 + 7.9109326a 3 
+ 0.8033946a 4 - [6.0758810 + 1.9532560a)JV/l a 2 
+ 0(a 3 ) (6.34) 

where clearly channel 1 correctly corresponds to the MOMq scheme definition as it is the only 
channel to contain the divergences in e. The quark external leg interchange also correctly 
emerges as a minor check on our computations corresponding to swapping p and q. From the 
renormalization constants in each scheme the numerical conversion functions from MOMq to 
MS for SU{3) are 

C™ 0Mq (a,a) = 1- [8.3578873 - 1.1720934a - 0.0820116a 2 - 0.5555556iV/] a 
- 131.2981279 + 7.8598968a - 0.3923795a 2 + 0.7219823a 3 

+ 0.0943409a 4 - [27.6257962 + 1.6277910a + 0. 1366860a 2 JiV/ 
+ 0.1543230^/1 a 2 + <3(a 3 ) 

C^ 0Mq (a, a) = 1 + [8.0833333 + 1.5000000a + 0.7500000a 2 - l.llllllliV/1 a 
+ [256.1034914 + 31.4182743a + 8.4375000a 2 + 2.8125000a 3 

- [53.9129277 - 1.6666667a] N f + 1.2345679 Nf 1 a 2 + 0(a 3 ) 
C^ 0Mq (a,a) = 1 + 3.0000000a 



r MOMq, 



+ 80.9357699 + 3.0725670a + 1.3465290a 2 - 5.9375000^ 
+ 0(a 3 ) 
a) = 1 — 1.3333333aa 



- 25.4642061 + 11.5753172a + 2.7222222a 2 - 2.3333333^ 
+ 0(a 3 ) . 

The analytic form of the first produces the coupling constant mapping which is 
«MOMq = °ms 



(6.35) 



+ 
+ 



52tt 2 - 78ip (!) + 993 C A - 240T F N f + 48^ '(§) - 32vr 2 - 432 

9768vrV( 1 r) - 7326(V>' (±)) 2 - 350910y/(±) + 2133^ w (i) 

+ 664848s 2 (f ) - 1329696s 2 (§) - 1108080s 3 (f ) + 886464s 3 (§) 

- 8944vr 4 + 2339407T 2 - 54432S - 164754C(3) + 3278628 



2 

_MS 
108 



+ 4617 



ln 2 (3)7r 



55404 



ln(3)vr 



4959- 



7T^ 



V3 v/3 VS_ 

+ \273m{4>' (\)f - 36480ttV(|) - 174816^'(±) + 1296<'(|) 

- 1959552s 2 (§) + 3919104s 2 (§) + 3265920s 3 (f ) - 2612736s 3 (§ ) 



+ 8704vr 4 + 116544^ 2 + 46656C(3) - 1960848 



13608 



ln 2 (3)7r 

~7T 



+ 163296 



In 3)vr t vr 3 
V + 14616^ 



CaCf 



24 



+ 



26112ttV(!) - 19584(V' / (i)) 2 + 295488^ 
- 373248s 2 (f ) + 746496s 2 (§) + 622080s 3 ( 



1,1 IV 

<3' 



7776^"' (|) 
-497664s 3 (§) 



+ 12032vr 4 - 196992vr 2 + 31104S + 435456C(3) + 256608 

3 



- 2592±^ + 31104^ + 2784^ 
v3 v3 v3 



+ [197568V>'(§) - 864^ w (i) + 186624s 2 (f ) - 373248s 2 (f ) - 311040s 3 (2 
+ 248832s 3 (§) + 2304vr 4 - 131712^ 2 - 217728C(3) - 1549440 



+ 1296 

+ U45447T 2 
+ ° (4s) 



ln 2 (3)7r 

~7T 



15552 



ln(3)7r 



1392 



7T 



C A T F N f + 115200T|A^ 2 



66816^(1) + 373248C(3) + 80352 C F T F N f 



3 

MS 
23328 



in the Landau gauge. Numerically we have 



(6.36) 



«MOMq 



+ 



16.7157746 - 2.3441868c* M g - 0.1640232a|jg - l.lllllllA^ 



"MS 



+ 



472.1590953 - 43.0575532a- 



0.7760123q4tct + 2.0207162a: 



+ 0.20885960^ - 



MS 



MS — MS 
83.1112168 - 0.6513959a 



MS 



MS 



+ 1. 2345679 Nf 



MS 



+ 



O (a- 



M 

MSJ 



(6.37) 



Thus, similar to the previous sections the explicit three loop Landau gauge renormalization 
group functions are 



/3 MOM >,0) = - 
+ 



YC A - 4 -T F N f 



a 2 - 



y C\ - 4C F T F Nf - ™C A T F N f 



[[71478(V/(i)) 2 - 953047r 2 ^(i) - 40266V'(i) - 7821V> w (i) 
- 2437776s 2 (f ) + 4875552s 2 (§) + 4062960s 3 (f ) - 3250368s 3 (f ) 
+ 52624tt 4 + 26844tt 2 + 199584S + 604098C(3) - 3593970 



16929 



ln 2 (3)vr 



+ 203148 



ln(3)vr 



7T 



+ 18183^ 



+ 



2069767rV(§) - 155232(V>'( ±)) 2 + 1951776^(|) - 4752<'(i) 

+ 7185024s 2 (§) - 14370048s 2 (§) - 11975040s 3 (f ) + 9580032s 3 (f ) 

- 563207T 4 - 1301184^ 2 - 171072C(3) - 159408 



+ 49896 



ln 2 (3)7r 

~7T 



- 598752 



ln(3)vr 



- 53592 



7T 

71 



c\c F 



+ 



34656ttV( 1 :) - 25992(^'(i)) 2 - 392328^' (§) + 6012^"' (|) 
+ 202176s 2 (f ) - 404352s 2 (§) - 336960s 3 (f ) + 269568s 3 (§) 
- 27584vr 4 + 2615527T 2 - 72576S + 578664C(3) + 3133080 

^3 



, dn 2 (3)vr t ln(3)vr 

+ 1404 y 16848- v ' 



V3 



+ 



J .is a- 



1508 



7T 



C z A T F N f 



88704(^'(i)) 2 - 118272 7 rV( 1 r) - 1387584^//(±) + 28512<'(£) 
+ 1368576s 2 (§) - 2737152s 2 (f ) - 2280960s 3 (§) + 1824768s 3 (f ) 



25 



- 366087T 4 + 9250567T 2 - 114048S - 1596672C(3) + 427680 
ln 2 (3K 



+ 9504- 



114048 



ln(3)7r 

7T 



10208^ 
>/3 



CaC f 



+ 



56448(^'(i)) 2 



75264tt 2 ( | ) - 592704?// ( § ) + 1 728^'" ( \ ) 
- 2612736s 2 (§) + 5225472s 2 (§) + 4354560s 3 (f ) - 3483648s 3 (f ) 
+ 204807T 4 + 395136vr 2 - 1306368C(3) + 1194912 



- 18144 



ln 2 (3)7T 

"7T 



+ 217728 



^ + 19488^ 



CaCfTpN^ 

+ [163584^/(1) - 1152<'(§) + 248832s 2 (£) - 497664s 2 (§) 



414720s 3 (f ) + 331776s 3 (§) + 3072vr 4 - 109056^ 2 
- 290304C(3) - 585792 + 1728 ^-^ 



- 20736 



ln(3)7r 

7T 



7T 



- 1856^ 
>/3 



C A T F N? 



+ 



430087rV( 1 r) - 32256(^'( 1 r)) 2 + 463104^' ( \\ 
- 497664s 2 (f) + 995328s 2 (f) + 829440^3(1, 
+ 13312vr 4 - 308736tt 2 + 41472S + 580608C(3) + 171072 
ln 2 (3)^ 



10368^ 
663552s 3 ( 



3456- 



+ 41472 



ln(3)vr 



+ 3712— 

v/3 



CpTpNf 



+ 



v/3 ' ™" v/3 
18432vr 2 - 27648^' (|) +497664C(3) - 352512] C F T F NfN { 



+ 0(a 5 ) 



23328 
(6.38) 



- [13CU - 8T F N f ] 



+ 



+ 



338vr 2 - 507V>'(|) + 864] C\ + [312^' (|) - 208tt 2 - 108] C A T F N f 
+ [312^(1) - 208^ 2 - 2808] C A C F 

2 

^128tt 2 - 192^(1) + 3024] C F T F N f ] ^ 

548808vrV(±) - 411606(^'( 1 r)) 2 + 4914^'(i) + 27729^ w (|) 
+ 8643024s 2 (f ) - 17286048s 2 (§) - 14405040s 3 (§) + 11524032s 3 (f ) 
- 2568807T 4 - 3276vr 2 - 707616S + 1016226C(3) - 2542266 



+ 



+ 6002li^ - 720252^ - 64467^ 
y/3 v/3 VS 



C\ 



+ 



+ 



745056^' (|)) 2 - 9934087r 2 ^'(i) - 8587296^(1) + 16848^'" (§) 

- 25474176s 2 (|) + 50948352s 2 (f) +42456960s 3 (f) 

- 33965568s 3 (|) + 286208tt 4 + 5724864tt 2 + 606528C(3) - 198288 

_ 176904^2^ + 2122848^^ + 190008^-11 C\C F 
\/3 v/3 v/3 

253296(V>' ( § )) 2 - 337728vrV( 1 r) + 1856304^ (§) - 28296^ 



in ( i 



2892672s 2 (|) + 5785344s 2 (§) + 4821120s 3 (f ) - 3856896s 3 (f ) 
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+ 1880327T 4 - 1237536tt 2 + 435456S + 3059856C(3) + 1166400 



20088 



I n 2 (3)tt 

"7T" 



+ 241056 



ln(3)7r 

7T 



. 7T° 



+ 21576^ 



C\T F N f 



+ 



>i> i i ■> 

3' 



499200vrV( 1 r) - 374400(V>' + 5997888^ (§) - 101088^'" ( 
- 4852224s 2 (f ) + 9704448s 2 (§) + 8087040s 3 (f ) - 6469632s 3 (f 
+ 1031687T 4 - 3998592vr 2 + 404352S + 5660928C(3) - 6368544 



- 33696 



ln 2 (3)vr 



+ 404352 



ln(3)vr 



7T° 



+ 36192— 



CaC% 



+ 



V3 ' ^3 

6113287rV(§) - 458496(V>'( 1 r)) 2 + 5660928^/(1) + 10368^'" (3) 
+ 15676416s 2 (§) - 31352832s 2 (f ) - 26127360s 3 (f ) 
+ 20901888s 3 (§) - 176128vr 4 - 3773952vr 2 

- 8584704C(3) + 7309440 + 108864 - - -y 

v3 



- 1306368 



ln(3)vr 

~7T 



7T*- 



- 116928^ 
>/3 



+ [6912<'(±) - 981504V' / ( 1 r) - 1492992s 2 (f ) + 2985984s 2 (f ; 
+ 2488320s 3 (f ) - 1990656s 3 (f ) - 18432vr 4 + 654336vr 2 

I 2 /q\ 

- 1244160C(3) + 155520 + 10368 - y 

v3 



+ 124416 



ln(3)vr 

7T 



7T 



+ 11136-= 

v/3 



C A T 2 F Nf 



+ 



230400ttV(±) - 307200(V' / (|)) 2 - 4188672^(1) + 62208V>"'( 1X 
+ 2985984s 2 (f ) - 597198s 2 (|) - 4976640s 3 (f ) + 3981312s 3 (§ ) 

- 63488vr 4 + 2792448vr 2 - 248832S - 3483648C(3) + 8118144 

.3 



ln 2 (3)vr ln(3)vr 
+ 20736 W 248832 



+ 



V3 ^/g 
165888^(1) - 1105927T 2 + 2985984C(3) 



22272-^ 



C F T F N f 



- 3608064] CfT^AT/ 



a 



139968 



+ 0(a 4 



(6.39) 



- \c A a 



+ 



+ 



26tt 2 - 39^' (|) + 90] C A + 36T F N f 



+ 



24V>'(|) - 16vr 2 - 216 



C F 



C A a 2 
72 



422167rV( 1 r) - 31662(^'(|)) 2 + 15066^(|) + 2133^"'( 1 ' 
+ 664848s 2 (f ) - 1329696s 2 (§) - 1108080s 3 (|) + 886464s 3 (f ] 
- 197607T 4 - 10044vr 2 - 54432S + 40338C(3) - 267786 



+ 4617 



ln 2 (3)7r 

~7T 



55404 



ln(3)7r 

~7T 



7T° 

4959^= 



+ 



57312(V''( 1 r)) 2 - 76416ttV(!) - 669600^/(1) + 1296^'" ( 
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- 1959552s 2 (f ) + 3919104s 2 (§ ) + 3265920s 3 (f ) 

- 2612736s 3 (f ) + 220167T 4 + 446400tt 2 + 46656C(3) + 66096 



13608 



In 2 (3)tt 

~7T 



+ 163296 



5^ + 146 16^ 
y/3 



C A C F 



+ 



38400(V>' (I)) 2 - 288007r 2 ^'(|) + 461376^) - 7776V>'"( i f ) 
- 373248s 2 (f ) + 74649652(1) + 622080s 3 (f ) - 497664s 3 (§) 
+ 7936vr 4 - 307584vr 2 + 31104S + 435456C(3) - 489888 



- 2592 



In 2 3)vr ln(3 vr m 1 vr 3 

\J— + 31104 V + 2784^ 

y/3 y/3 V3 



cl 



1036807$ AT/ 



+ [145152^(|) - 864^ w (i) + 186624s 2 (f ) - 373248s 2 (f ) 

- 311040s 3 (f ) + 248832s 3 (f ) + 2304vr 4 - 96768tt 2 - 15552C(3) 



+ 273456 + 1296 



ln 2 (3)7T 
V3 



- 15552 



+ 



230407T 2 - 34560^' (§) + 264384] C F T F N f 



V3 y/S 



C A a 3 
31104 



C A T F N f 



+ 0(a 4 ) (6.40) 



and 



7^ 



MOMq 



(a, 0) = [25Ca - QC F - 8T F N f ] 



C F a 2 



+ 



r 2 



+ 
+ 
+ 

+ 



3900^(|) - 2600vr 2 - 13230C(3) + 10287 

2224vr 2 - 3336^' (|) + 5184C(3) + 160921 C A C F 



832vr 2 - 1248-0 (!) + 3456C(3) - 8784 C A T F N f 
768^(1) - 512vr 2 - 7776] C F T F N f + 1152T| AT/ 

C F a 2 



576V» (|) - 384vr^ - 4536 



C% 



432 



+ 0(a 4 



(6.41) 



The corresponding arbitrary linear covariant SU (3) expressions are 

/3 MOMq (a,a) = - [11.0000000 - 0.6666667N/ ]a 2 

- 102.0000000 + 15.2372141a - 1.3839787a 2 - 0.4920696a 3 

12.6666667 + 1.5627912a + 0.2186976a 2 ] N f 
1843.6527285 + 422.0731852a + 123.3734958a 2 - 19.5130255a 3 

- 3.5055190a 4 - 0.0961314a 5 

588.6548455 + 60.5454812a + 16.3955703a 2 + 0.9282359a 3 

- 0.0000059a 4 ] N f + 22.58781 18 A/] a 4 + 0(a 5 ) 

[0.6666667AT/ - 6.5000000 + 1.5000000a]o 

- 46.6391320 + 22.5608759a - 6.2001294a 2 + 0.8789652a 3 



MOMq/ \ 
1 A (a, a) 



9.4117058 + 1.5627912a - 0.3906512a 
2027.7437143 + 333.3082218a + 184.2382915a 
+ 12.6718858a 4 + 1.9200786a 5 



24.3519716a d 
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- 415.6990152 + 49.4053081a + 9.7003844a 2 - 3.7908552a 3 



0.4783683a 4 



+ [11.1788079 - 1.3021707a] Nf 



0(a 4 



7 



MOMq 



(a, a) 



[0.7500000a - 2.2500000]a 

- 13.2020072 + 12.3112512a - 2.5140879a 2 - 0.6855174a 3 
- 0.7500000TV/] a 2 

- 740.1341645 + 1.8665778a + 100.6450352a 2 - 3.4355918a 3 



8.7678000a 4 - 1.0965129a 5 



N f + 2.50000007V f 2 



MOMq/ v 
7^, M (a,a) 



75.5032720 + 4.1186469a + 1.7109768a 
+ 0(a 4 ) 
1.3333333aa 

22.3333333 + 2.4900784a + 8.1255824a 2 + 1.2186976a 3 

2 



'/ 



+ 



+ 



- 1.33333337V/] a 2 
341.8989103 + 182.9132891a + 43.9801061a 2 + 74.9283461a 3 
+ 21.4352687a 4 + 1.9493562a 5 



- 52.1916907 - 3.1076278a - 2.1669462a 2 
+ 0(a 4 ) . 



TV 



f 



0.88888897V 



a' 
(6.42) 



Unlike the previous two sections we close this section with a few remarks on other renormal- 
ization group functions in QCD. Recently, there has been interest in the RI'/SMOM renormal- 
ization scheme which was developed in |49j . Briefly the scheme was introduced to renormalize 
flavour non-singlet quark currents in a quark 2-point function at the symmetric subtraction 
point. The motivation for such a scheme is that it circumvents potential infrared singularities 
of the measurement of the same Green's function on the lattice in the chiral limit for excep- 
tional momentum configurations, [20] . More specifically the scalar, vector and tensor currents 
have been considered at one and two loops in this scheme, |49} 1431 150] . and the associated 
conversion functions for each operator relative to the MS scheme have been determined. Thus 
the RI'/SMOM three loop operator anomalous dimensions have been deduced in an arbitrary 
covariant gauge. More recently, the full set of amplitudes has been provided for each of these 
currents at two loops in [51] . Equally the same information has been provided for the n = 2 
and 3 moments of the flavour non-singlet Wilson operator central to deep inelastic scattering, 
\52\ 153] . As an extension of the RI'/SMOM analysis we have determined the three loop anoma- 
lous dimensions of the scalar and tensor quark currents as well as the n = 2 moment of the 
Wilson operator in the MOMq scheme. Numerically the SU(3) values are 

7 M° Mq (a, a ) _ _ 4.0000000a - [7.5709424 + 8.3450254a + 0.3121855a 2 + 1.79187637V) 

+ [324.9490278 + 84.8639648a - 37.0694357a 2 - 1.1713670a 3 



0.2647488a 4 - [l6.2843867 + 18.7759399a + 0.5878184a 2 



N 



2.66666677V, 



a 3 + 0(a 4 



(6.43) 



MOMq/ \ 
7 T H (a, a) 



1.3333333a + 



20.3014252 + 8.1573041a + 1.8959382a 2 - 0.58789317V/ 
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+ 



276.1622997 + 149.4337760a + 83.3978821a 2 + 32.6152122a 3 



+ 3.1873834a 4 
+ 2.07407417V/ 



61.3516342 + 4.4008934a + 0.1928562a" 



a 3 + 0(a 4 ) 



N f 



(6.44) 



and 

MOMq/ \ 



3.5555556a + 



45.2685593 + 15.7105153a + 3.0417366a 2 - 2.62643457V, 



+ 



818.7944256 + 154.0549414a + 59.1199636a 2 + 6.2844353a 3 



+ 0.3234930a 4 



+ 4. 73923797V, 



140.9726462 + 5.3687957a + 0.8615923a' 



a 3 + 0(a 4 



N f 



MOMq/ \ 
Tw 2 ,12 ( a > a ) 



1.7777778a 



16.5859286 + 7.5010250a + 0.9582634a 2 - 0.7715739/V/ 



298.1271156 + 48.6705419a + 27.5758422a 2 + 3.1422177a 3 
+ 0.1617465a 4 - [54.9281680 + 0.1449714a + 0.2531120a 2 



N f 



+ 2.17522117V 



/ 



a 3 + 0(a 4 ) 



MOMq/ \ 



0(a 4 



(6.45) 



where S and T respectively denote scalar and tensor currents. The Wilson operator second 
moment three loop anomalous dimension corresponds to the 11 entry of the dimension two 
upper triangular mixing matrix, [52]. The 22 entry corresponds to the total derivative operator 
for this operator moment but its anomalous dimension is equivalent to that of the vector current, 
[52j . As the vector current is conserved and hence a physical operator its anomalous dimension is 
zero to all orders in perturbation theory consistent with the Slavnov- Taylor identity. Whilst we 
have included a subset of the operators considered in recent RI'/SMOM analyses for comparison 
we have not recorded the complete set for the MOMq or either of the other MOM schemes 
discussed here. The reason for this is simple. In deriving (|6.43p . (|6.44p and (|6.45p we have 
arrived at the results by constructing the respective conversion functions and then using the 
analogous relation to (|3.9p for operator renormalization. It transpires, like our observation for 
the wave function and gauge parameter conversion functions, that the respective conversion 
functions are formally the same as those already given in [491 H3"l [50] for the quark currents and 
that for W2 given in |52j . The only difference here in producing (|6,43p . (|6.44p and (I6.45P is that 
the coupling constant mapping between the schemes is not trivial as it is in the RI'/SMOM case. 
Therefore, if one is interested in the structure of the three loop anomalous dimensions in the 
MOMggg and MOMh schemes it is merely a simple exercise to construct them from the known 
conversion functions. We should add that in deriving the two loop parts of (I6.43p . (I6.44p and 
(|6,45p we have not only used the renormalization group construction but also carried out the 
explicit two loop renormalization of each of the three operators in the MOMq scheme. Therefore, 
the two loop terms of each have been derived directly and also indirectly from the respective 
one loop conversion functions which we have constructed explicitly. 



7 Discussion. 



We close with some remarks and give an overall perspective on our computations. First, we have 
computed the two loop structure of each of the trivalent QCD vertices with a linear covariant 
gauge fixing at the symmetric subtraction point for the MS scheme as well as the three MOM 
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schemes. Although we have only given the MOMggg, MOMh and MOMq versions for each of 
the respective triple gluon, ghost-gluon and quark-gluon vertices, the values of say the triple 
gluon vertex in the MOMq scheme, if of interest, can be readily deduced from the results given 
here via the Slavnov- Taylor identities of QCD, |21j . In addition we have deduced the basic 
three loop renormalization group functions including the /3-functions in each of the three MOM 
schemes in an arbitrary linear covariant gauge. As there has been numerical estimates for some 
of these quantities in the Landau gauge in [22] we need to comment on how our results compare 
in addition to earlier remarks. First, we consider the SU(3) Landau gauge two loop mapping 
of the coupling constant from the MOMi schemes to MS. Specifically, the polynomials in Nf 
of the perturbative map were given in [22] and the coefficients were denoted by dy where 
/ is the loop order and j is the power of N f in the polynomial. As d^ om , d^° Mi and d^ om 
can all be computed analytically without approximation and which we agree with for MOMggg 
and MOMq, we focus on the remaining two, d^ om and d^ 0Mi , which could only be evaluated 
numerically in |22j . In Table 3 we give the results of |22] in its notation, in order to preserve the 
error estimates, together with the numerical evaluation of our analytic results for SU(3) in the 
same convention. Aside from the MOMh scheme the agreement is remarkably close. However, 
that for (fM° Mh differs by about 5% whilst that for d^ OMh is significantly different. We believe 
that this is a transcription inconsistency in the presentation of the results in [22J akin to that 
noted earlier here and in [28j for the amplitudes of this vertex. This is also because we do not 
tally with the one loop coefficients (i^o and ci^fOMh which are known exactly, [21] , We have 
checked that we obtain the same one loop MOMh ghost-gluon vertex counterterm given in [21] 
exactly for all a. Equally we have obtained the same vertex counterterm as [2T] for the quark- 
gluon vertex as a check on our overall procedures which were the same for each vertex. For 
instance, the one loop master integral denoted by / in |21] and [22] arises in 

d MOMh It derives 

directly from the vertex counterterm only and the coefficient of / is therefore in a one-to-one 
mapping to the coefficient of / in diff^^. Comparing the counterterms of our results and [22], it 
is clear that the coefficients differ but we find our expressions are consistent with [21]. Therefore, 
the one loop agreement of our MOMh results with [21] suggests that we have a more credible 
MOMh coupling constant mapping as we have followed the same symbolic algebraic algorithm 
as that for MOMggg and MOMq. 





,MUMggg 
"20 


,MUMggg 
"21 


^MOMh 


^MOMh 


JVlOMq 
"20 


jJVLOMq 


Ref [22] 
This paper 


59.8(8) 
60.02892 


- 12.6(2) 
- 12.63031 


35.88(2) 
40.12125 


- 5.0707(6) 
- 5.34725 


29.53(1) 
29.50994 


- 5.1961(4) 
- 5.19445 



Table 3. Comparison of two of the two loop SU(3) coupling constant relation coefficients for 
each scheme in the notation and convention of |22j . 



Concerning the final Landau gauge St/ (3) /3-functions of [22] . we have already briefly noted 
the accuracy of our exact expressions in [21] . Summarizing, for the coefficient of the Nf indepen- 
dent three loop term the agreement was around 2% for the MOMggg and MOMh /3-functions 
and less than 0.2% for MOMq. In Table 4 we have presented the coefficients of the polyno- 
mial in A^- of the three loop term of the /3-function in each scheme. These are denoted by fiij 
where again I is the loop order and j is the power of Nf in the polynomial. Only the MOMggg 
scheme has a cubic term in A^ and it is known analytically. So the precise numerical agreement 
merely reflects this. The linear term in Nf in the MOMggg scheme is virtually zero because 
of an accidental cancellation for the colour group SU(3) which was noted in [22]. This is the 
reason why a large error was recorded for this coefficient though we obtain the same sign for 
our coefficient. Overall for MOMggg and MOMq all the coefficients are very close to the central 
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values given in [22]. That for /3 30 had a relatively large error, possibly as a result of the 

purely gluonic graphs which might have a more slow convergence in the asymmetric expansion 
parameter approach used in |22j . However, it is still reassuring that our result emerges so close 
to the central value in this particular case which roughly has a 2% error. Therefore, in this con- 
text the slightly smaller error on /3^ OMh suggests that ultimately the MOMh /3-function of [22] 
has been correctly obtained. Although the agreement of the other Nf polynomial coefficients 
is reasonably decent there does not appear to be as close an overlap compared to say those 
for MOMq. This is probably, because, as we noted earlier, the amplitude for the triple gluon 
and quark-gluon vertices for the channel containing the poles in e, or equivalently the channel 
corresponding to the vertex Feynman rule, are in very good agreement with the Landau gauge 
estimates given in |22j . Therefore, overall we are confident that our results are correct. 



Scheme 


Article 


/3so 


031 


/3 32 


033 


MOMggg 


Ref [22] 
This paper 


24(2) 
24.5466309 


0.04(63) 
0.0088426 


- 1.05(3) 
- 1.0482740 


0.0415330 
0.0415330 


MOMh 


Ref [22] 
This paper 


44.82(5) 
43.9608273 


- 9.730(5) 
- 9.6507368 


0.3276(1) 
0.3359815 




MOMq 


Ref [22J 
This paper 


28.86(3) 
28.8070739 


- 9.206(3) 
- 9.1977320 


0.35322(7) 
0.3529346 





Table 4. Comparison of coefficients of three loop Landau gauge /3-function in the convention of 

[22]. 

One final point worth emphasising in this regard is to note that when the computations of [22] 
were being carried out the technology available in terms of algorithms, such as [3D], and computer 
power was not comparable to that of current levels. Therefore it is a remarkable achievement 
that overall the results are in solid argreement. This is partly because the expansion technique of 
|22j to evaluate a 3-point Green's function at the symmetric point in terms of 2-point functions 
to allow the application of Mincer, requires a large amount of integration by parts especially 
as the powers of the propagators increases in each term of the expansion. Despite this, however, 
it would seem that the method of [22] could be extended in principle to tackle the construction 
of the vertex functions to three loops. Whilst there are now three loop anomalous dimensions 
and /3-functions available as intermediate checks on any numerical estimates, it would provide a 
reasonable estimate for four loop renormalization group functions including the /3-functions in 
physical MOM schemes. Such higher order computations may be necessary soon in order to not 
only assist with lattice matching of vertex functions but also for more accurate determination 
of the value of the strong coupling constant. 



Acknowledgement. The author thanks Prof. K.G. Chetyrkin for discussion on |22] and 
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A Tensor basis and projectors. 

In this appendix we list in turn the basis of tensors used as well as the coefficients for the 
projection of each of the three vertices. Away from the symmetric subtraction point the tensor 
basis is enlarged. 
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A.l Triple gluon vertex. 



For the triple gluon vertex we have the basis tensors 



>ggg 

(4)/ii/o 



' ! 1/1/ 



T)lxvV<J 



V: 



(p, q) 



V 



ggg / 
(10) fiua- ^ 



\P»P»P* , 'Pf^AP'l) = J^PuPvQ" > ^(P-?) = ^iP^vPa 



i 1 



(13)nva 



(p,q) 



V, 



(11)\we 



(p,q) 



1 



2%Pv1° 

(A.l) 



The first six tensors represent the structures which appear in the triple gluon vertex Feynman 
rule. As there are fourteen tensors for this vertex then in order to simplify the presentation we 
have partitioned the projection matrix into a 3 x 3 block matrix with partitions of sizes 6, 4 
and 4. Defining 



A^ggg 



MW 



27(d - 2) 

then in d-dimensions the sub-matrices for the projection are 



A4ff s 
M ggg 



M 



33 



/ 36 




M ss. 



11 



V 
/ 





36 


18 







36 



18 



M 



13 



12 
12 

48 
24 
24 
24 



48 
12 
12 
24 
24 
24 



M 



22 



M 



23 



M 
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18 



12 
48 
12 
24 
24 
\ 24 

( 64(d + l) 
32(d + 1) 
32(d + l) 
V 32(d + l) 

/ 16(d + 4) 
8(4d + 1) 
8(4d + l) 
8(d + 4) 

/ 12 48 12 

12 12 48 

48 12 12 

\ 24 24 24 



V 



18 

18 



36 

36 



24 \ 
24 
24 
48 
48 

48 J 

32(d + 1) 
32(2d - i; 
16(d + 1) 
16(d + l) 

16(d + 4) 
8(4d + l) 
8(d + 4) 
8(4d + 1) 

24 24 

24 24 

24 24 

48 48 



\ 



18 





36 / 



M ggg 



M ggg 



1 48 


24 


24 


24 \ 


48 


24 


24 


24 


48 


24 


24 


24 


24 


48 


12 


12 


24 


12 


12 


48 


V 24 


12 


48 


12 J 


48 


24 


24 


24 \ 


24 


48 


12 


12 


24 


12 


12 


48 


24 


12 


48 


12 J 



32(d + l) 32(d + l) 

16(d + l) 16(d + l) 

32(2d-l) 16(d + l) 

16(d + l) 32(2d-l) 

16(d + 4) 8(d + 10) \ 

8(d + 4) 16(d + 4) 

8(4d + l) 16(d + 4) 

8(4d + l) 16(d + 4) J 

24 \ 
24 
24 

48 J 



) 



33 



M 



32 



M 
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in d-dimensions. 



( 16(d + 4) 8(4d + l) 8(4d + l) 8(d + 4) \ 

16(d + 4) 8(4d + l) 8(d + 4) 8(4d + l) 

16(d + 4) 8(d + 4) 8(4d + l) 8(4d + 1) 

V8(d+10) 16(d + 4) 16(d + 4) 16(d + 4) ^ 

/32(2d-l) 16(d + l) 16(d + l) 32(d + l)\ 
16(d + l) 32(2d-l) 16(d+l) 32(d + l) 
16(d+l) 16(d + l) 32(2d-l) 32(d + l) 
\ 32(d+l) 32(d+l) 32(d + l) 64(d + l) j 



(A.2) 



A. 2 Ghost-gluon vertex. 

At the symmetric subtraction point there are two basis tensors for the ghost-gluon vertex which 
are 

■)CCB 



ccg 
(2)0 



1c 



Hence, in this case the projection matrix is relatively simple and is given by 

A4 ccg = 

in d-dimensions. 



1 

3 



4 2 
2 4 



(A.3) 



(A.4) 



A.3 Quark-gluon vertex. 

The quark-gluon vertex tensor basis involves six independent tensors. These are the same as 
those used for the vector current insertion in a quark 2-point function in [51] and are given by 



la 



-r>qqg i \ PuP nqqg I \ P^n 



2 ( 3 ) 



(A.5) 



where the generalized d-dimensional 7-matrices denoted by r|*\'"' Xn were defined earlier. Equally 



■2 > V V>^ 
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16(d- 1) 
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1) 


8(d- 


1) 


4(d + 2) 





1 


6 


8(d- 1) 


4(4d- 


-7) 


4(d- 


1) 


8(d- 1) 





36(d - 2) 


6 


8(d- 1) 


4(d- 


1) 


4(4d- 


-7) 


8(d- 1) 







12 


4(d + 2) 


8(d- 


1) 


8(d- 


1) 


16(d- 1) 







V 


















-12 J 



the same projection tensor in d-dimensions emerges as in |51| which is 



M qqg 



The partition of the projection matrix into the and r^? 7 sectors is apparent. 
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